High order correlation functions for self interacting scalar field in de Sitter space 
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^ I (Dated: Notes du 30 Octobre 2003) 

^ , We present the expressions of the three- and four-point correlation functions of a self interacting 

■ light scalar field in a de Sitter spacetime at tree order respectively for a cubic and a quartic potential. 

^ ' Exact expressions are derived and their limiting behaviour on super-horizon scales are presented. 

Q ' Their essential features are shown to be similar to those obtained in a classical approach. 

00 ' I. INTRODUCTION 

T— I , 

^-H ■ The computation of high order correlation functions of scalar fields in a de Sitter spacetime can be of interest 
^ ' for investigations of the physics of the early universe. There are indeed a growing number of indications that in 
, its early phase the universe underwent an inflationary period '3| that can be accurately described by a de Sitter 
phase during which the energy density of the universe is thought to be dominated by the self energy of a scalar 
' field Tj. Its fluctuations, or the fluctuations of any other light scalar fields, are thought to be the progenitors of the 
large-scale structure of the universe 3] . The statistical properties of the induced metric fiuctuations then depend on 
I the potential landscape in which the fields evolve. It is unlikely that the fluctuations produced along the inflaton 
. direction can be significantly non-Gaussian as stressed in recent works 0, 0- It is possible however that if there 
exist other self interacting light scalar fields during that period primordial non Gaussian metric fluctuations can be 
generated |E S B HI- The details of the effects induced by such self- interacting fields is based on the statistical 
1 properties of the metric fluctuations that quantum generated field fiuctuations can induced. 
O ' The purpose of this paper is to investigate the foundation of these calculations by calculating the correlation 

' properties of a test scalar field in a de Sitter background. 
^ [ The content of the paper is the following. In the second part we present the basis of such computations. In the third 
• • 1; we explore the expression of the leading order term of high order correlators for different self-interacting potentials. 
^ ' In particular we give explicit results for the superhorizon limit that corresponds to modes that can be observed today. 
Finally we give insights on the physical interpretations of those results. 



II. CORRELATION FUNCTION COMPUTATIONS 



A. Free field behavior 



For the purpose of our calculations we assume that the infiationary phase can be described by a de Sitter background 
epoch. This description is only approximate since during the inflationary period the Hubble constant slowly varies 
with time but it is a framework in which all calculations can be pursued analytically. We have checked that to a large 
extent the results would not be affected if the background is changed Q . 

A de Sitter spacetime in flat spatial section slicing is described by the metric, 

"^'^ " {ikiy + -JydxMx^") (1) 

which is conformal to half of the Minkowski spacetime. The conformal time is related to the cosmological time by 



(2) 
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and runs from — oo to 0, the limit r] representing the "infinite future". 

For a minimally coupled free quantum field of mass m, the solution can be decomposed in plane waves as 

vo{^,7j) = J d^k [vo{k,fi)bue'''-^ + v*o{k,rj)ble''''-^] (3) 

where we have introduced v = ax, a hat referring to a quantum operator. In this Heisenberg picture, the field 
has become a time-dependent operator expanded in terms of time-independent creation and annihilation operators 
satisfying the usual commutation relations [6k, ^j^.,] = (5Dirac(k — k'). We can then define the free vacuum state by the 
requirement 

6k |0) = for all k. (4) 

As it is standard while working in curved spacetime the definition of the vacuum state suffers from some 

arbitrariness since it depends on the choice of the set of modes xo(fcj?7)- They satisfy the evolution equation 

^ - Z2 — No = 0, (5) 



with |c2p — |cip = 1, where H^^^ and 



the general solution of which is given by yfwr]J\ ciHi}\kri) + C2H^\k'q) 

are the Hankel functions of first and second kind and with = 9/4 — rri^/H^. Among this family of solutions, it 
is natural to choose the one enjoying the de Sitter symmetry and the same short distance behavior than in Minkowski 
spacetime. This leads to 

vo{Kv) = ^^H^^\k7j). (6) 

This uniquely defines a de Sitter invariant vacuum state referred to as the Bunch-Davies state vacuum In the 
massless limit, the solution (O reduces to 

This result gives the expression of the equal-time two-point correlator of the Fourier modes, 

(X(ki)x(k2)) = (5Dh-ac(kl+k2)P2(fcl) (8) 
7:r2 2 / -i \ 

P2(fcl) = -J- l + _ . (9) 



2fci V kfr]^ 

A remarkable result is that in the superhorizon limit {kri — > 0) is that the phase of vo{k,r]) freezes and it reads, 

i;o(fc, 77) = --!=. (10) 

Note that the Fourier modes of x are frozen only in the case of a massless field in a de Sitter background. As a result 
the field x behaves like a classic stochastic field with fluctuations whose 2-point correlator is, 

(X(ki)x(k2)) = ^ <5Dirac(kl + k2) (11) 

And since x is a free field, its superhorizon fluctuations follow a Gaussian statistics. 



B. Computation of higher-order correlation functions 



Having determined the free field solutions, one can then express perturbatively the iV-point correlation functions 
of the interacting field, X: in terms of those of the free scalar field. The equal time correlators we are interested in 
are expectation values of product of field operators for the current time vacuum state. Such computations can be 
performed following general principles of quantum field calculations [ij. The simplest formulation it is to apply the 
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evolution operator U(r]Q,r]) backward in time to transform the interacting field vacuum into the free field vacuum at 
an arbitrarily early time tjq so that, 

(i-ki.-.^^kj = m-\7jo,v) vi.i---Vk„ C/(%,^)|0) (12) 

where |0) is here the free field vacuum |15|. It is implicitly assumed in this expression that the coupling of the field 
X is switched on at time tiq. We will see in the following that the choice of 770 is not important as long as it is much 
earlier than any other times intervening in the problem. 

The evolution operator U can be written in terms of the interaction Hamiltonian, Hj, as 



(13) 



If one is interested only in a single vertex interaction quantity, the evolution operator can be expanded to linear 
order in Hj, 

U{m,v) = Id-i r dij' Hj{rj') (14) 

so that the connected part of the above ensemble average at a time rj finally reads, 

(«ki...^^kjc--i /'d77' (0|[^;kl...Wk„,i^/(^?')]|0) (15) 
•J no 

where the brackets stand for the commutator. 

The result is expressed in terms of the Green function 

G(fc, ry, rj') ^ -L (^1 _ (^1 + exp[ifc(,7' - 77)] (16) 

defined as (0|?;(k, 77)z;(k', 77')|0) (5Dirac(k + k')G(/c, 77, 77'). 

Before we proceed to explore explicit cases, let us note that as long as calculations are restricted to tree order, the 
very same calculation can be done assuming that x is a classic stochastic field whose stochastic properties are initially 
those of the quantum free field. Not surprisingly one finds the same formal expressions! 

III. EXACT RESULTS 

Results can be given in a closed form for simple self-interacting potentials. We give in this section explicit results 
for cubic and quartic potentials. Because of its renormalization properties, quartic potential is a natural choice to 
consider. We will however see that the case of a cubic potential can be relevant when finite volume effects are taken 
into account (see Ref. Ts^ for more details on finite volume effects). 

In the following we thus assume that Hi is of the form 



(17) 



and we are then interested in the computation of the leading order part of the connected part of the ensemble average 
of products of p Fourier modes of the fields. 
In case of a quartic coupling, we have, 

(i-ki ■■■vu,)c = -iA5Dh.ac(ki + . . . + k4) y" d77' [G{ki ,v,v')--- G{h, V, v') - G* {h , 77, 77') • . . G* (^4, 77, 7;')] (18) 

This expression is unfortunately quite cumbersome to compute. It depends, for symmetry reasons, on the norms of 
the four wave vectors fci,...,fc4 in the following combinations 

^1 = ^fc. (19) 

i 

Y^hk, (20) 

hkjkk (21) 



7r2 



i<j<k 

7r4 = "Y 

i<j<k<l 
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shape dependence of f 
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FIG. 1: Dependence of the function as a function of ki for fixed values of fci, k2 and fcs. The solid line corresponds to 
configuration fci=0, k2 = \, ^3=1, ^4; the dashed line to = ^2=1, ^3 = ^, ^4 and the dotted line to ki — ^, ^^2 = ^, ^3 = ^, ^4- 



The expression (jl8|l can then be rewritten as 



TTl — 2 TTl 7r2 — TTl TTs + 3 7r4 + ( — TTl TTs + 3 TTl 7r2 TTs — TTl 7r4 — 3 7r2 11^) 7] + 3 7r4 

-TTi (tti^ - 3 TTl 7r2 + 3 tts) I [(l - + J?'^) cos(7ri Tj) + T] (tti - TTa T]'^) sin(7ri ry)] ci(-7ri 77) 

- [r] (tti - TTs r]'^) cos(7ri 77) + (-1 + 7r2 77^ - 714 ry'^) sin(7ri 77)] si(-7ri 77)| j (23) 

where the Coslntegral (ci) and Sinlntegral (si) functions are defined, for 77 < 0, by 

dn' 

exp[— 1/077'] = ci(— fc77) + isi{—kr]). (24) 

-00 

The superhorizon limit (i.e. kii] ^ I for all i = 1 ... 4) of Eq. H23|l is 

{vki---VUi)c = '^ff'T ^'-^ liTi'^ - 2 7ri^7r2 - tti 713 + 3 714 - tti (7ri^-3 7ri 7r2 + 3 7r3) [7 + log(-7ri 7;)] |(25) 
24 7r4 7ri 77* L ^ ' J 

where 7 is the Euler's constant (7 « 0.577). This implies that the 4-point correlator of the actual fields x reads 

In this expression, terms of the order of kii] have been neglected. This result illustrates the transition to the stochastic 
limit. The function C is an homogeneous function of the Fourier wave- numbers 

/■ni ^\ -TTl'' + 2 7ri^ 7r2 + 7ri 7r3 - 3 7r4 

C[{h}) = 5 —: ^ • (27) 

TTl (7ri-' — 3 7ri 7r2 + 3 7r3) 

The dependence of the function ( on the wavelength ratios is illustrated on Fig. ^ It is found to be relatively 
weak. We did not fully explore its properties but one can explicitly show for instance that if one of the wave vectors 
vanishes, C is smaller than —2 and reaches its minimum ,—8/3, for a symmetric configuration of the three remaining 
wave vectors. From the plot on Fig. ^ it appears that the minimum value of ^ is —3.2 and is reached for a "square" 
configuration, that is when the four wavelengths are of equal length. The maximum value, —2, is reached when two 
of the wavelengths vanish. The overall variations of C with the wavelength ratios are therefore rather mild. It is then 
legitimate to describe the four-point correlation function as the sum of products of two-point power spectra, e.g., 

(Xki . . . Xkjc = fchac (J2 ^4({fcJ) (28) 

mh}) = ^3({fcJ)5]nii- (29) 
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square configuration 
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FIG. 2: Behaviour of the function Q4 as of function of time. The transition to the superhorizon behavior (dashed line) is 
shown. The function Q4 is shown here for a "square" configuration (fci — k2 = = k4) as a function of ^477 = ^ kirj. 

elongated rectangle configuration 
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FIG. 3: Same as Fig. |21for a "rectangular" configuration {ki — k2 — 4:, k-^ = 4^4). 




The vertex value, is given by 

^3({fc.}) = ^ [7 + Cih) + log (30) 

that carries a weak geometrical dependence with the wave vectors geometry through ^. 

When the term log {—rjj^ ^i) is large (and negative), that is when the number of e-folds, N^, between the time of 
horizon crossing for the modes we are interested in and the end of inflation is large, the vertex value is simply given 

which corresponds exactly to the value that was obtained in Ref. ^] from the classical evolution of the stochastic field 
on superhorizon scales (see comments in last section). 

The result obtained above is however more complete since it contains next to leading order terms. It allows, for 
instance, to estimate the validity regime of the superhorizon result. It shows in particular that the mode coupling 
induced on subhorizon scales are negligible as soon as the number of e-folds after horizon crossing exceeds a few units. 
This is made clear when one considers the reduced four-point correlation function, Q^, defined as, 

^ ,fu ^ \ - P4{ki,k2,k3,ki) 

P2(fci)P2(fc2)P2(fc3) + sym. ^^^^ 

where P2 is defined in Eq. (Q. The behaviour of (54 as a function of time and for different configurations of the 
wavevectors is depicted on Figs. (I2l4fl and one can convince hismself that it converges rapidly toward the superhorizon 
result as soon as — logJ^kiV ~ ^(1)- 

Similar results can be obtained for the 3-point function in the case of a cubic potential, ^(x) — In this 

case the quantity to compute is 

r -dn' 

(^^ki-.-i^ka) =-iA<5Dirac(ki + ... + k3) J [G(A:i , 7?, ??') • ■ • ^(fca , 77, " G* (fci , ry, 7?') . • . G* (fcg , r?')] • (32) 
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equilateral triangle configuration 




which gives, 



{Vki . . . Wk: 



Log[ki ri] 

FIG. 4: Same as Fig.|5|for a "triangular" configuration (fcl — 0, k2 — ks — ^4) 



■ Suhnc (ki 



12 TT^ if H 

TTi^ - 2 TTi 7r2 - TTa + ("TTi^ TTg + 3 TTa TTg) 7?^ 

- (tti^ - 3 TTi 7r2 + 3 tts) I [(1-71-2 ?7^) cos(7ri 77) + 77 (tti - tts 77^) sin(7ri 77)] ci(-7ri 77) 

- [77 (tti - TTs 77^) cos(7ri 77) + (-1 + TT2 V^) sin(7ri 77)] si(-7ri 77) |^ (33) 

To a factor 277, one recovers basicahy the same expression as for the fourth cumulant H23|) when is set to zero. In 
the superhorizon limit, the expression of this cumulant reads 



(Xki ■•■Xk3> = - 



where (^3 is simply given by 



XH^ '^Dirac (Ek.) 

12 Ukf 

C3({fc»}) =C(fcl,fe,fc3,fc4-^0) 



To express it in another way we still have 

(Xki ■■•Xks) 



'^Dirac(Ek.) P3i{k^}) 



p,{{h}) = -3({fc.})Eni3 



(34) 

(35) 

(36) 
(37) 



with a vertex value 7^3 to be taken in the appropriate limit, fc4 — s- 0. The shape dependence of the vertex, or the time 
dependance of this cumulant thus reproduces that of the fourth cumulant as shown on Figs. ^ (solid line) and^ 



IV. COMMENTS 



We have obtained some closed forms for the 3- and 4-point correlation functions of a test scalar field with a 
self-interacting potential of order respectively 3 and 4 in a de Sitter background. 

It is interesting to note, as it was already in the literature 's^, that the superhorizon behaviour of the field can be 
obtained from a simplified Klein-Gordon equation for the field evolution on superhorizon scales 



dV 

X + 3Hx = -^(x) 
ax 



(38) 



solved perurbatively at first order in A. In this equation x has to be understood as the filtered value of x a fixed 
scale that leaves the horizon at a given time to- The previous equation is then valid for t > only where the field 
can be described by a classical stochastic field. Moreover in writing this equation one also makes the assumption that 
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its r.h.s. can be computed from the filtered value of the field (which is not necessarily identical to what would have 
been obtained from a filtering of the source term). The equation H38|) can be solved perturbatively in A. At zeroth 
order x is constant and, at first order, it reads 

if x^iP is the leading order value of the field at horizon crossing. Note that t — to can be rewritten as N^/H where Ne 
is the number of e-folds since horizon crossing. It implies that if the number of e-folds is large enough the term x'o'^ 
should become negligible. As a consequence the leading order expression of the first nontrivial high order cumulant of 
the one-point PDF of x takes either the form 3 i^six^)'^ or 4 1^3 (x^)^ for respectively a cubic or a quartic potential. The 
value of the coefficient that enters these expressions is precisely the one obtained in Eq. H30|l in the superhorizon 
limit, that is 1^3 = —XNf./{3H^). And indeed if one had to compute such cumulants from either expression (|29|) or 
(|37|l the integration over the wave vectors would have led to the very same expressions in the superhorizon limit. That 
shows that the late time behavior of the cumulants we found comes in fact from a simpler dynamical evolution. It 
actually demonstrates that the computations in the classical approach sketched here can be put on a firm ground. It 
also allows investigations of more subtle effects, such as the finite volume effects where more than one filtering scales 
have to be taken into account, that cannot be properly addressed in a classical approach. 

Regarding finite volume effect it is interesting to investigate how the results for the cubic and the quartic potentials 
could be related together. We refer here to our companion paper "iS"! where we give a comprehensive presentation 
of finite volume effects on observable quantities. Suffice is to say that even if one assumes x has a self interacting 
quartic potential, a non-zero third order correlator might be observable. The reason is that one cannot have access 
to genuine ensemble averages but to constrained ensemble averages such as (xki • ■ •Xka)^ which is the expectation 
value of Xki ■ ■ ■ Xka for a given value of %, average value of x over the largest scale available in the survey in which 
the correlators are computed. 

To leading order in x with respect to its variance, a-^, this constrained average value reads 

X 

(Xki • ■ -Xk,,)- = {Xvki ■ ■ -Xka x)c (40) 

If, for i — 1...3, the wavelengths 1/ki are much smaller than the survey size then the expression of 
(Xki • • • Xk3 x)c x/f^-l is that of (xki ■ • ■ Xka) when the self-interacting field x evolves in the potential Axx'^/3!, which 
is the cubic term in x in the expansion of A (x + x)^/4! • This shows that the two cases described in this paper are 
consistent with one another and that the functional relation between the two had to be expected. 

To conclude we have established in this paper a number of generic results that put the tree-order computation of 
higher order correlations in a de Sitter background on a secure ground. We leave for further studies the examination 
of more realistic cases in which the background expansion is more complex. 
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